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We present a quenched lattice QCD calculation of spin-1/2 five-quark states with uudds quark 

^- 

content for both positive and negative parities. We do not observe any bound pentaquark state in 

<— | 

these channels for either I = or I = 1. The states we found are consistent with KN scattering 

Oh. states which are checked to exhibit the expected volume dependence of the spectral weight. The 

results are based on overlap-fermion propagators on two lattices, 12 3 x 28 and 16 3 x 28, with the 

>: 

same lattice spacing of 0.2 fm, and pion mass as low as ~ 180 MeV. 
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I. INTRODUCTION 

Since the reported discovery [l] two years ago of an exotic 5-quark resonance, named 
Q + (uudds), with a mass of about 1540 MeV and a narrow width of less than 20 MeV, there 
has been a rapid growth of interest in the subject. Eleven more experiments have reported 
the observation of the state J2J. It also stimulated the search for other pentaquarks J3|. 
It should be pointed out, however, that there are also ten experiments reporting negative 
results 4|. One has to wait for high statistics experiments to clarify the situation in order 
to establish the exotic state beyond doubt. 

The strangeness quantum number of 9 + is S = +1, but its isospin and spin-parity 
assignments are undetermined by the experiments. Based just on the valence quark content, 
the isospin could be 0, 1, or 2. The spin-parity could be ^ , | , or higher. The isospin 
would have to be established by discovering the other charge states, while the spin-parity 
assignment will have to await detailed measurements of decay angular distributions. 

The experiments were inspired by the Skyrme model prediction [5J], and the experimental 
discoveries have in turn spawned intense interest on the theoretical side, with studies ranging 



from chiral soliton and large N c models |6[, quark mode 



p, quark models [7] 



J9(, and recent lattice calculations |l0j. 



QCD sum rules |9(, and recent lattice calculations 
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KN phase-shift analysis 



II. INTERPOLATING FIELDS 

Unlike ordinary mesons (qq) and baryons {qqq), pentaquarks do not have a unique color 
structure aside from being a color singlet. For a spin-1/2 pentaquark state of the type uudds, 
there are two decay modes with different isospin content, K°p and K + n. The simplest local 
interpolating field can be written as a color-singlet configuration of a product of color-neutral 
meson and baryon interpolation fields, 

X T = e a b c (uT* Cld b ) [U C {S e lb d e ) T{U~ d}] , (1) 

where sum over all the color indices {a, b, c, e} is implied. The minus sign is for isospin 1=0 
and plus sign for 1=1 respectively. A slight variation with a different color contraction is 
given by 

x t = e a*c (u Ta C l5 d b ) [u e (s e l5 d c ) T {u ~ d}\ , (2) 



where the color indices e and c are positioned differently. Both interpolation fields in Eq. ([!} 
and Eq. (J2J) have been used in a lattice calculation to study the pentaquark uudds [lOj. 
Another possible 1 = interpolation field inspired by the diquark-diquark-antiquark picture 

is 

X l = e 9ce e gfh e abc (^Ta^^ (yTf Crd h\ p^-ljTe 

= e abc (u Ta C l5 d b ) (u Tc CTd e ) TC- 1 ^ 

_ e abc ^Ta^^ (V 6 CTcf ) FC^S^, (3) 

where V = {S,A} = {1,7^75}, and we used the antisymmetric tensor relation e 9ce e 9 ^ h = 
5 c ^5 eh — 5 ch 5 e f to obtain the second part of the equation. The interpolation field in Eq. (121) 

n n 

with T = S has been studied in Ref. yjj and Ref. |12l |. 

These interpolation fields couple to states with both parities. For X3, as for the standard 
nucleon interpolation field, the positive-parity state propagates in the forward (backward) 
time direction in the upper (lower) Dirac components of the correlation function, while the 
negative-parity state propagates backward (forward) in the upper (lower) components. On 
the other hand, for Xi and X2, the positive-parity state propagates in the forward (backward) 
time direction in the lower (upper) Dirac component of the correlation function, while the 
negative-parity state propagates backward (forward) in the lower (upper) component. In 
our calculations, we use both upper and lower components to improve the statistics. 

In general, since X\iX2 and x\ have the same quantum numbers (e.g. 1 = and J p = 
1/2 1 * 1 ) they should project out the same states, albeit with different spectral weights. It 
is interesting to note that they can be explicitly related. Indeed, despite their apparent 
different color-spin structures, the KN type interpolation fields and that of the diquark- 
diquark-antiquark type are related by a factor of 75 and a Fierz re-arrangement which 
switches the roles of the u quark and s in Eq. (JTJ, leading to the following expression 

(u Ta C l5 d b ) l5 u c (s e l5 d e ) = -(u^C^d^iu^Cd^C- 1 ^ 

+h K u Ta C 1 ,d b ){u Tc C 1 ^) 1 ,C- l s Te 
-Uu^C^d^u^Ca^d^a^C- 1 ^ 

o 

+\(u Ta C 1 sd b )(u Tc C 1 , l5 d e ) 1 , l5 C- 1 s Te 
+^u Ta C 75 d b )(u Tc C l5 d e ) l5 C- 1 s Te . (4) 



In this expression, a sum over dummy indices /i and v is implied, so the right-hand-side 
has 16 terms. We use the 7-matrix relation {7^,71/} = 26^. We see that the first term 
on the right hand side is just the first term in the second part of Eq. (J3J). A similar Fierz 
transform of X2 m Eq. (J2J) will lead to the second term in the second part of Eq. (jHJ). Since 
75 multiplication reverses parity, the extra 75 factor in front of the u quark on the left hand 
side serves to match the explicit parity of both sides. Therefore, schematically one can write, 

75 x (KN interpolation field) 
= \ x (diquark-diquark-antiquark interpolation field) 

+ other terms, (5) 

where the KN interpolation field is either xi or X2- More precisely, 

7s (xT - X2) =^(x! + Xi) • (6) 

The fact that the KN interpolation field and the diquark-diquark-antiquark interpolation 
field are directly related further enhances the argument that both interpolation fields couple 
to the same physical spectrum with different strengths. This exercise also suggests other 
possibilities for the diquark-diquark-antiquark interpolation fields for pentaquarks. In view 
of this, it is not surprising that the ground state results using KN type interpolation fields |lOl | 
largely agree with those |ll| using the diquark-diquark-antiquark type. But it is a puzzle 
that Ref. |l2|, using the same diquark-diquark-antiquark type interpolation field, pro duces 



airk type 1: 
and [ll| 



Although Ref. 



12j uses 



qualitatively different results than those of Ref. 

the overlap fermion, which has better chiral properties than the Wilson fermion adopted 

in Refs. H| and |ljj, one would not expect a qualitative difference in these calculations 



between the two fermion formalisms for the pion masses larger than 440 MeV. We will 
compare our results with these lattice calculations later. 

The correlation function is obtained by Wick-contractions of all possible quark pairs. In 
the case of the interpolating field in Eq. (P), it has four terms. Due to isospin symmetry in 
the u and d quarks, the two diagonal terms are equal, and so are the two cross terms. The 
zero-momentum correlation function without the upper(lower) Dirac component projection 
reads 

G 5q (t) = £< X(x) X(0) > = £( X(x) X(0) )diag ± ( X(x) X(0) ) cross , (7) 



where the diagonal contribution is given by, 



Y,(x(x)x(0))dL 9 = 2Y, e *bce a >b>c> x 



S% (u)Tt 



S(d)tf(s)] tr [S bb '{d)S cc '{u) 



S a a Uu)tr V bb (d } s } d)S cc (u) 



+ 



S aa ^gt* (Qgcc ^ 



ad 



Tr 



S(d)S\i 



S aa (u)V bb (d } s,d)S cc (u) 



a/3 



2 2_j € abc^a'b'c' x 



and the cross contribution is 

E(x(x)x(0))cLs 

V^(d,s,u)tr[S bb '(d)S cc \u) 
+ [S aa '(u)S bb '(d)V cc '(d,s,u) 
+ [v aa '(d,s,u)S bb '(d)S cc '(u) 

- ~S aa '(d)V bb '(d,s,u)S cc '(u) 



ap 
ap 



(8) 



(9) 



In the above expression, S = S q (x, 0) is the fully-interacting quark propagator, S. = 
(CSC" 1 ) 7 with C = C75 and V(u,s,d) = S(u)S^(s)S(d). The trace on spin-color is de- 
noted by 'Tr' while 'tr' represents the trace on spin only. In Eq. (J7J) plus sign is for isospin 
1 = and minus sign for I = 1 respectively. 

III. RESULTS AND DISCUSSION 



Our results below are obtained on two lattices, 12 3 x 28, and 16 3 x 28, using the Iwasaki 

Ql — I 
and the overlap fermion action [14(. The lattice spacing of a = 0.200(3) 

fm was determined from /„- |l9| for both lattices, so the box size is L = 2.4 fm and L = 
3.2 fm, respectively. Our quenched quark propagators cover a wide range of quark masses, 
corresponding to pion mass from 1293(20) MeV down to 182(8) MeV. The strange quark 
mass is set by the meson, and corresponds to a pseudoscalar mass r% ~ 760 MeV. In Fig. 
n we plot the nucleon and the kaon masses as a function of rri^. A naive linear extrapolation 
to the physical pion mass yields a kaon mass which is within ~ 7% of the corresponding 
experimental value (star symbol). The Iwasaki gauge action is an 0(a 2 ) renormalizat ion- 
group improved action which allows the use of relatively coarse lattices without suffering 
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FIG. 1: Nucleoli and kaon masses as a function of m^ for 3.2 fm lattice. Experimental values are 
represented by the star symbols. 



from large discretization errors. The overlap fermion action preserves exact chiral symmetry 
on the lattice, and thus it has no 0(a) error. One further finds that 0(a 2 ) errors are 
small for the meson masses lal and renormalization constants llfil ll7L 181. Owing to the 

nn 

relatively gentle critical slowing down [la, \\j\ , it allows us to work at unprecedented small 
quark masses. The relatively large box size ensures that the finite-volume errors are under 
control. At our lowest pion mass, the finite- volume error is estimated to be ~ 2.7% [l9(. We 
have used the combination of Iwasaki gauge action and overlap fermion action with local- 
local correlators in a number of recent studies, including chiral logs 19] and baryon excited 
states [201 . To handle the excited states, we use the recently developed constrained-curve- 
fitting algorithm - the sequential empirical Bayes method 
here contain 80 gauge configurations each. 



2111. The two ensembles studied 



A. KN scattering states 

It is worth noting that the entire five-quark spectrum is contained in the correlation func- 
tion in Eq. (JJJ), including a tower of KN scattering states and possible bound pentaquarks. 
Considering the presence of a pentaquark state, the above correlation function can be written 
as a sum of exponentials 



G 5q (t) 
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(10) 
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FIG. 2: Correlation functions for the positive parity (first column) and the negative parity (second 
column) for 2.4 fm (first row) and 3.2 fm (second row), respectively (at a pion mass around 200 
MeV). One should notice that the scale in the top left figure is different than the others. This 
shows that for the positive parity channel the ghost state is more prominent for the smaller lattice, 
while no ghost state is detected for the negative parity channel. 

It should be stressed that the ordering of possible pentaquark states and KN scattering 
states with energies E l KN in Eq. (|l()jl is not known a priori, and must be determined by 
fitting the data. Furthermore, one needs to discern the nature of the fitted states in order to 
distinguish if they are KN scattering states or bound pentaquark states. The parameter E l KN 
denotes the KN two-particle energies with zero total momentum. Since the KN interaction 
is relatively weak, their values are expected to be near the two-particle threshold energy 
denned as 

E K (p) + E N (p) = ^m 2 K +p 2 K + ^Jm 2 N +p 2 N . (11) 

We use the discrete momentum available on the lattice: PK{ n ) — Vn(2/a) sin(7r/L) for the 
kaon and PN{n) = ^/n{l/a) sin(27r/L) for the nucleon, where n = 0, 1, 2, • • •. As an example, 



TABLE I: The KN threshold energies corresponding to the first few discrete momenta at three 
different pion masses are given for the two lattices. In each block, the first column is for the 
12 3 x 28 lattice, the second column 16 3 x 28 lattice. 



n 


E KN (GeV) 


E KN (GeV) 


E KN (GeV) 




(nm = 


182 MeV) 


(m* = 


438 MeV) 


(m n = 


692 MeV) 





1.47 


1.47 


1.77 


1.77 


2.16 


2.16 


1 


1.80 


1.67 


2.06 


1.94 


2.41 


2.30 


2 


2.07 


1.85 


2.30 


2.10 


2.63 


2.44 


3 


2.31 


2.00 


2.52 


2.24 


2.83 


2.56 



the KN threshold energies corresponding to the first few discrete momenta are given in 
Table HI On our lattice, the energies calculated from the above forms of lattice momenta 
deviate about 2-3% from the corresponding energies calculated from Pk = Pn = 2n/L. 

We will refer to these discrete states as p — 0, p — 1 for n — 0, n — 1, and so on, keeping 
in mind that their actual values are given in Table [I] One can see the expected down-shift 
of the states on the larger volume. The spacing between neighboring states decreases as n 
increases and the momentum states are somewhat more packed for heavy quark cases. The 
discrete KN scattering states play different roles in the positive-parity and negative-parity 
channels. In the positive-parity channel, they are in a relative P-wave, so the spectrum 
starts at p = 1, which has a raised threshold of 1.80 GeV at m^ = 182 MeV, while in the 
negative-parity channel, they are in S-wave and the spectrum starts at zero momentum 
with a threshold of 1.47 GeV. If there exists a pentaquark near 1.54 GeV, it would lie below 
(above) the KN threshold in the positive (negative) parity channel. This means that it 
would be much easier to extract it in the positive-parity channel than in the negative-parity 
channel. 

The experimental values for the KN scattering length (volume) are given in Table ILT1 
The numbers are taken from Ref. 22j. The interaction is very weak in the 1=0 channel 
(with a small attraction in the P-wave), and is slightly repulsive in the 1=1 channel. The 
5-wave scattering length on the lattice is related to the energy shift in a finite box of length 



TABLE II: Experimental values for the KN scattering length (volume). 



scattering 



S-wave (fm) 
P-wave (fm 3 ) 



1=0 1=1 



0.0 ±0.03 -0.32 ±0.02 
0.08 ±0.01 -0.16 ±0.1 



Lby 



pint / p , p \ _ 27ra o 



l±c 1T + c 2 - 



+ 0(P" 6 ), (12) 



HknL 3 

where a is the S-wave scattering length and hkn is the reduced mass of the KN system. 
Using the experimental numbers and neglecting c\ and C2, the estimated energy shift for 
I = 1 is about 8 MeV on our 16 3 x 28 lattice, and about 18 MeV on the 12 3 x 28 lattice. 
As it will become clear later, these values are consistent with our results. One could turn 
this argument around and use our simulation results to extract information on the KN 



scattering 



24J . a subject outside the focus of this work. 



B. Ghost states 

Since we are working in the light quark region with the quenched approximation, there 
is an additional complication arising from hairpin diagrams corresponding to the would- 
be 7]' meson which is degenerate in mass with the pion. The resulting states containing 
the would-be rf are called ghost states since they are unphysical quenched artifacts with 
negative spectral weights. Because of the negative weight the key signature for the presence 
of ghost states is the negative correlation function at very small pion mass. Such quenched 
artifacts have been observed in both the meson [25] and baryon J20J sectors. In the case 
of pentaquarks, the ghost state is NKr] 1 which has intrinsic positive-parity with relative S- 
waves among the hadrons. In the l/2~ channel, there is a relative P-wave between a pair of 
N, K, and the ghost rf. The threshold of this state is thus raised on the lattice: for example, 



Eknt)' = Ek + En + E v > = Jm 2 K + p 2 K + m^ + Jm 2 + v\ where p starts from n = 1 for a 
P-wave between K and //. On our smaller lattice, the value is about 2.23 GeV at 188 MeV 
pion mass. So the ghost state in the l/2~ channel lies quite high while the lowest states are 
expected to be the KN scattering state in S-wave threshold and the possible pentaquark. 
In the l/2 + channel, the situation is the opposite. The K, N, and the ghost 7/ are all in 
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FIG. 3: Correlation functions for the positive parity (first column) and the negative parity (second 

column) for 2.4 fm (first three rows) and 3.2 fm (second three rows) lattices, respectively. These 

correlation functions are for three pion masses (in MeV) . One should notice that there are negative 

dips in the positive parity channel which indicate the presence of ghost states. These KNr]' ghost 

states decouple from the correlation function at pion mass above ~ 300 MeV. 
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relative S'-wave which has a threshold of EKNrfin = 0) = mK+mN+m n . Using experimental 
values for pion, nucleon and kaon Eknti'{ti = 0) = 1.57 GeV, which is right near the possible 
pentaquark of 1.54 GeV. So in the l/2 + channel, the ghost state lies relatively low and plays 
a significant role, just like in the S -n (1 535) channel where the ghost state lies lower than the 



Su when the quark mass is light 20] (i.e. with pion mass lower than 300 MeV). So in this 
channel, the ghost state and the potential pentaquark state are candidates for the ground 
state, followed by the excited state of n — 1 KN scattering state. In our algorithm, the 



ghost state is modeled as —W g (l + -E^e - ^ 4 where E n = Jm% + pi and E g is constrained 
near its threshold value J20J. The fact that the ghost state has a negative spectral weight 
—W g is crucial for its isolation from the non-ghost states. 

In Fig. |2l we present 1 = 1 correlation functions (at a pion mass around 200 MeV) for 
both parity channels and for two lattice volumes (2.4 and 3.2 fm). The top two figures are 
for 2.4 fm and bottom two are for 3.2 fm. A row-wise comparison of these figures reveals that 
the positive parity correlation functions have a negative dip which indicates the presence 
of quenched ghost states. On the other hand, correlation functions for the negative parity 
channel are always positive which confirms the absence of low-lying ghost states in that 
channel. Comparing the left two figures (one also has to notice the difference of scales) it is 
clear that the ghost states in the smaller volume are more prominent. 

In Fig. 01 we plot a few more 1 = 1 correlation functions. As in Fig. |21 the left side figures 
are for positive parity and the top six figures are for the smaller volume lattice. Figures 
on the left side show the effect of ghost states in the correlation function. Note that the 
ghost states are only noticeable in the very light quark mass region. The effect of the ghost 
states decreases as the pion mass increases, and decouples from the correlation function near 
pion mass around 300 MeV. The right side figures are for the ne gati ve parity channel and 
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12|, the pion masses 



they are always positive. In the previous lattice calculations [lOj, 

are above 440 MeV, and therefore, there is no need to consider the ghost states. Since our 

lowest pion mass is 182 MeV, we will need to take the ghost states into account in the l/2 + 

channel. 
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FIG. 4: The computed ground state mass in the I{J P ) = 0(l/2~) channel as a function of m^ 
for the two lattices L=2.4 fm and L=3.2 fm. The curve is the KN threshold energy in the S-wave 
ExN( n = 0) = rriK + mjv- The bottom figure is an enlarged version of the top figure for the small 
quark mass region. 

C. Negative-parity channel 

In Figs. |U and EJ we show the results for the ground state mass as a function of m^ in 
the I(J P ) = (1/2 - ) and 1 (l/2~) channels respectively. Also plotted is the KN threshold 
energy in the S'-wave Ekn^ = 0) = tjik + itln which is the same on both lattices. There 
is no need to consider ghost states in these channels, which is supported by the fact that 
the correlation functions are positive throughout. The ground state energy on the smaller 
lattice (L=2.4 fm) is consistently higher than that on the larger one (L=3.2 fm). At the 
lowest mass, the energy coincides with the S"-wave threshold, meaning that there is little 
interaction, which is consistent with the experimental fact of zero scattering length (see 
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FIG. 5: The computed ground state mass in the I(J P ) = 1 (l/2~) channel as a function of m^. 
As above, the curve is the KN threshold energy in the S'-wave EKN{n = 0) = rriK + niN, and the 
bottom figure is an enlarged version of the top figure for the small quark mass region. 

Table HU). 

Here we point out that the uncertainty in the strange mass determination (which we set 
by the meson and resulting in our kaon mass being higher than the experimental value 
by about 7%) is harmless. We note that changing the kaon mass affects both the data 
and the energy threshold simultaneously, so that the difference between them is essentially 
unaffected. 

We notice that at larger m^, the fitted ground state is higher than the mK+rriN threshold. 
In the higher quark mass region, the excited states have larger weight relative to the ground 
state. This, together with a somewhat smaller mass gap between different momentum states 
at higher quark masses, results in an effective mass which continues to drop at the largest 
time slices. As a result, in the higher quark mass region it is more difficult to isolate the 
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function of m^ (for 3.2 fm lattice). The first excited states (square) are shown along with the 
ground state (circle). The first few scattering states corresponding to different discrete momenta 
(for n = 0,1,2,3 and 4) are shown by dashed lines. 

ground state, and the fitted lowest state is a mixture of the ground state and the excited 
states with n = 1 and higher. The fact that the fitted lowest state energy on the smaller 
lattice (L=2.4 fm) is higher than that on the larger one (L=3.2 fm) for the same higher mj is 
consistent with this interpretation since the lattice momenta on the smaller lattice is larger 
than that on the larger lattice, leading to a mixed state with higher energy in the smaller 
lattice. In the low quark mass region, the excited states have smaller relative weights, and 
thus we are able to fit the ground state more accurately close to the mass threshold (i.e. 
m,K + rax with n — 0). However, the lack of statistics prevented us from resolving the 



n 



1,2 and 3 states which are closely packed together 



2l|. As can be seen from Table |TJ 



the separation of KN states with different momenta ranges from 150 to 200 MeV which are 
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FIG. 7: The computed mass in the I{J P ) = 0(l/2 + ) channel as a function of m^ for the two 
lattices L=2.4 fm and L=3.2 fm. The two lower curves are the KN threshold energies in the 
P-wave Ekn{ u = !)• The two higher curves represent the energies of the non-interacting ghost 
states. 

much lower than the usual separation between the ground states and their radial excitations 
in ordinary hadrons. The latter are usually about 500 - 600 MeV. Consequently, our fitted 
first excited states appear to be a mixture of the n = 1 and higher momenta states. 

As far as the ground states are concerned, our results more or less agree with those of 

Y2\. It is noted in Ref. |ipj 



Ref. llOj and ll|, but disagree qualitatively with those of Ref. 
and 11] that they have seen a low- lying excited state above the KN mass threshold and 
they interpret it as the pentaquark state. From the top figure of Fig. we see that the 
first excited state near the chiral limit is ~ 1.85 GeV which is substantially higher than 
1.54 GeV, the experimentally observed pentaquark state. As explained above, we interpret 
it as the mixture of the n — 1, n — 2, and possibly n = 3 KN scattering states. We 
observe similar behavior for the 1=1 case also (bottom figure). There is no candidate for 
a pentaquark between the KN threshold (n = 0) and the first KN scattering state (n = 1) 
which is contrary to the finding of Ref. [ll| . We tried to accommodate an extra low-lying 
pentaquark state in between our ground and first excited states. However, the x 2 fit always 
rejects such an intermediate state. This implies that our data do not favor a pentaquark 
state in between the two scattering states with lowest momenta. This, however, does not 
preclude a pentaquark state nearly degenerate (i.e. within 100 MeV) with the KN scattering 
state at threshold as was observed in Ref. [lTj with a two-channel approach. We shall study 
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FIG. 8: The computed mass in the I(J P ) = 1 (l/2 + ) channel as a function of m^ for the two 
lattices L=2.4 fm and L=3.2 fm. The two lower curves are the KN threshold energies in the 
P-wave Ekn{ u = 1)- The two higher curves represent the energies of the non-interacting ghost 
states. The bottom figure is an enlarged version of the top figure in the lower quark mass region. 

this possibility in the future |26J . 

D. Positive-parity channel 

The fitted ground state mass as a function of m\ in the I(J P ) = (l/2 + ) and 1 (l/2 + ) 
channels are shown in Fig. Eland Fig. |H1 respectively. In these channels, the ghost NKr]' must 
be included, as discussed above. In the fitting model, the NKt]' ghost state, pentaquark, 
and KN P-wave scattering state are considered. We found a ghost state and a KN scattering 
state, but not a pentaquark state near 1.54 GeV. In the lower energy 1 = 1 channel, we have 
tried to see if our data could accommodate three states, but the x 2 /dof would simply reject 
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it. The energy of the KN scattering state lies higher on the smaller lattice (L=2.4 fm) than 
that on the larger lattice (L = 3.2 fm). This mainly reflects the fact that the momentum 
corresponding to n = 1 is larger on the L = 2.4 fm lattice than that of on the L = 3.2 
fm lattice. At the lowest mass, the energies almost coincide with the P-wave thresholds, 
meaning that the KN interaction is weak, consistent with experiment. At higher m^, we 
again observe that the ground state lies substantially higher than the P-wave KN threshold 
energy. We believe this is a result of a mixture with several higher momenta states, as 
proposed in the last section to explain the similar behavior for the negative parity states. 
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FIG. 9: A comparison of ground state masses as a function of rn^ for 1 = and 1 = 1 channels 
both for the negative (top) and the positive (bottom) parities. 

In Fig. El we compare the ground state energies for 1=0 and 1=1 on the 3.2 fm lattice. 
We observe that in the \ channel, although within errorbars, the energy of the 1=1 state 
tends to be higher for all masses. This tendency is inverted for the ~ channel. Both 
these observations are consistent with Ref. [lfj. Our conclusion that the pentaquark state is 
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absent below the KN P-wave threshold again agrees with those of Ref . |l0|, LL3J and disagrees 
with that of Ref. Il 



E. Volume dependence 

In a box, lattice states have a volume dependence of Jl/V (where V is the spatial 
volume) from the normalization on the lattice. For a one-particle state, a point-source and 
zero momentum point-sink correlation function is given by 



Git) = 5>lx(a?,*)*(0)|0> 



l<0|x(0)|n)P 



■e 



-M n t 



^ 2M n V 

Y,W n e- Mnt (13) 



where 



H-«* (14) 

2M n y ' 

is the spectral weight corresponding to the mass M n , which has no explicit volume depen- 
dence. On the other hand, for a non-interacting two-particle scattering state, with a total 
zero momentum, the corresponding correlation function can be written as 

G(t) = £(0|xi(x,£) X2 (x,£)xi(0)x2(0)|0> 

= 1 y l(0|Xi(0)|ni>| 2 |(0|X2(0)|n 2 >| 2 ,,-^+^t 

mn 2 2E ni V2E n2 V 

= <y W ni W n2 c _ {Eni+Enr))t ^ 

nw2 

which has an explicit inverse volume factor. Though the total momentum between the 
two particles is zero they can have finite lattice momentum, and thus, in Eq. E3 we used 
energies E ni and E n2 instead of masses M ni and M n2 . Since the KN interaction is very 
weak, approximating the spectral weight of the KN scattering state with the non-interacting 
expression (Eq.ll5j) should be a reasonable one, as far as the volume dependence is concerned. 

Besides the explicit volume dependence, the spectral weight in the noninteracting case i.e. 

w„,w„ 



/ ni '► no 



v has implicit volume dependence from W ni and W n2 . However, in our case, the 
spectral weights for the nucleon, Roper, and Sn have very small volume dependence (at 

18 



a few percent level) between the 12 3 x 28 and 16 3 x 28 lattices [2Jj. As a result, we can 
expect that even for a worst case the combined volume factor from W ni and W n2 would be 
much smaller than the explicit volume factor V, which is 2.37 for our lattices. Therefore, 
it can be concluded that the spectral weight W for a one-particle state should not exhibit 
explicit volume dependence whereas a two-particle scattering state will show an explicit 
volume dependence in its spectral weight. 
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FIG. 10: Ratio of spectral weights on the two lattices in the negative parity 1 = (top figure) and 
1=1 (bottom figure) channels. The line at 2.37 is the expected volume dependence of the spectral 
weight. 

This volume dependency of spectral weight can be exploited to check whether the state 
extracted is a one-particle bound pentaquark state or a two-particle KN scattering state. 
This is especially important to check for the negative parity channel where the threshold 
scattering states are close to the presumed pentaquark state. On the two lattices we used, 
the ratio of the spectral weights is expected to be W(12)/W(16) ~ 16 3 /12 3 = 2.37. Fig. [TU] 
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FIG. 11: Ratio of spectral weights of the two lattices for I(J P ) = 1 (l/2 + ) channel. The line at 
2.37 is the expected volume dependence of the spectral weight. 

shows the ratio for the negative channel. Indeed it deviates significantly from unity, and 
the value is consistent with ~ 2.37 within error bars. This is strong supporting evidence 
that the observed state is a KN scattering state. The situation is similar in the positive 
parity channel, as shown in Fig. ^2 We conclude from the volume dependence that the 
ground states we observe (besides the ghost states) are the KN scattering states, not the 
pentaquark states. Similar volume studies have been carried out earlier 20j to verify that 
the Roper and Sn states and the ghost Nr/' states do have the expected one-particle and 
two-particle volume dependences respectively. Therefore, we caution that before concluding 
that a state near the KN scattering threshold is a pentaquark state, one must study the 
volume dependence of its spectral weight. If the weight has a large volume dependence one 
should not mistakenly identify it as the physical pentaquark state. The presence of KN 
scattering states complicates the study of possible pentaquarks in the 5-quark spectrum. It 
is therefore desirable to have a correlation function which is completely dominated by the 
two-particle scattering states and compare it with the point-source and point-sink correlators 
which have both the two-particle scattering states and the possible one-particle pentaquark 
state. One can also project out the ground state and excited state explicitly by using multi- 
operator cross correlators (similar to Ref. |ll|). The details of such analysis will be presented 
elsewhere J26|. 
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IV. CONCLUSION 

Our results based on the overlap fermion with pion mass as low as ~ 180 MeV seem 
to reveal no evidence for a pentaquark state of the type uudds with the quantum numbers 
I(J P ) = (0,1) (| ) near a mass of 1540 MeV. Instead, the correlation functions are domi- 
nated by KN scattering states and the ghost KNr]' states in the l/2 + channel at low quark 
mass (pion mass less than 300 MeV). Our results are consistent with the known features of 



the KN scattering phase-shifts analysis [22|. We have checked that the observed KN states 
exhibit the expected volume dependence in the spectral weight for two particles in a box. 
We advocate the use of this volume dependence to uncover the character of the states found 
in multi-quark calculations on the lattice. 

Our conclusion is in contradiction with the other lattice calculations which have claimed 



a pentaquark signal of either negative parity 10, 11], or positive parity [12J, in the vicinity 
of 1.54 GeV. 
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